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Abstract 

The  purpose  of  this  report  is  to  examine  the 
concepts  upon  which  diffusion  theory  is  based  and,  in 
particular,  to  show  that  the  rate  of  diffusion  of  a 
conducting  fluid  across  a  magnetic  field  is  indepen- 
dent of  the  field  If  the  proper  scaling  is  used. 
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Introduction 

The  problem  of  the  diffusion  of  an  ionized  fluid 
across  a  magnetic  field  has  been  investigated  both  the- 
oretically and  experimentally    ,   The  aim  has  been  to 
determine  how  the  coefficient  of  conduction  or  diffu- 
sion is  affected  by  the  presence  of  a  magnetic  field 
B.   Both  theory  and  measurements  have  shown  that  the 
rate  of  diffusion  perpendicular  to   B   is  decreased, 
this  being  attributed  to  a  decrease  in  the  coefficient 
of  conductivity.   The  latter  is  usually  defined  as  the 
coefficient  of  proportionality  between  the  current   J 
and  the  electric  field  E.   This,  as  we  shall  show  later* 
is  an  ambiguous  statement  and  the  dependence  of  conduc- 
tivity on  B   depends  on  how  exactly  conductivity  Is  de- 
fined.  The  question  as  to  which  inverse  power  of  B   the 
conductivity  depends  on  experimentally  has  been  settled 

TTT 


Chapman  and  Cowling.   The  Mathematical  Theory  of 
non-uniform  gases.   Cambridge  University  Press,  1939 

A.  Guthrie  and  R.  K.  Wakerllng.   The  Characteristics 
of  Electrical  Discharges  in  Magnetic  Fields.   (McGraw- 
Hill  Book  Co.,  Inc.,  New  York,  1949)   National  Nu- 
clear Energy  Series,  Flutoniura  Project  Record,  Vol,  5 
Division  1. 

A.  Simon.  Ambipolar  Diffusion  in  a  Magnetic  Field. 
Physical  Review,  Vol.  98,  317-318,  1955. 
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1      (2.) 
only  recently  by  Simon  in  favor  of  a  —2  law    • 

B 
In  this  paper,  the  problem  will  be  treated  from 

the  view  point  of  continuum  theory,  most  existing 
theories  of  diffusion  having  been  based  on  kinetic 
theory.   The  purpose  is  to  elucidate  some  of  the  vm- 
derlying  concepts  rather  than  to  present  explicit  for- 
mulas for  diffusion  rates.   These  will  be  given  in  a 
note  to  follow  since  they  might  be  of  some  use  in  ex- 
perimental work.   In  particular,  we  shall  show  that 
when  the  equations  of  the  fluid  are  expressed  in  di- 
mensionless  form,  only  one  essential  parameter  ^ 
remains.   If  a  strict  comparison  is  made  between  ex- 
periments, this  parameter  should  be  kept  constant. 
We  shall  show  that  as  a  consequence,  the  velocity  of 
diffusion  is  independent  of  B.   If,  however,   P  is 
allowed  to  vary,  i.e.,  strict  scaling  is  no  longer 
maintained,  then  a  dependence  on  B  can  be  obtained. 
For  one  special  type  of  scaling,  the  velocity  of  dif- 
fusion will  be  shown  to  be  proportional  to  — p  for 

B"^ 
large  B,   in  agreement  with  experiment. 

The  conclusion  that  the  coefficient  of  conduc- 
tivity is  independent  of  B  is,  at  first  sight,  rather 
surprising.   However,  it  does  not  contradict  any  pre- 
viously established  theories  or  experimental  results. 

Rather,  it  serves  to  emphasize  that  this  coefficient 

__ —  ■         - 

^^Cf.  reference   (1), 
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has  no  precise  meaning,  independent  of  the  definition 
which  is  assigned  to  it  or  the  experiment  which  ineasiires 
it.   The  same  situation  exists  in  regard  to  the  notion 
of  "mean  free  path"  which  may  have  many  different  mean- 
ings depending  on  how  it  is  measured.   In  order  to  il- 
lustrate this  point:   Ohm's  law  may  be  written  as 
j»  =  yE'   +  6E'xB  where   y  and  6   depend  on  B,  J* 
and  E'   being  the  current  and  electric  field  as  viewed 
from  a  coordinate  system  moving  with  the  velocity  of 
the  fluid.   If   Y  is  identified  with  the  coefficient  of 
conductivity,  we  may  state  that  the  latter  depends  on  B. 
However,  if  the  formula  is  inverted  with  E'   expressed 
in  terras  of   J',   then  the  coefficient  of   J'   turns  out 
to  be  independent  of  B.   Another  example  which  could  be 

used  arises  from  the  equation  of  motion  for  B,   namely, 

2 
^  +  (u-V)  B  +  B(Vu)  =  ^V  B 

where   u  Is  the  velocity  of  the  fluid,   ^t,   the  per- 
meability, and  4,   the  conductivity  in  the  absence 
of  B.   Apart  from  the  convection  terms,  the  equation 
has  the  form  of  a  diffusion  equation  for  B  -  hence,  we 
state  that  in  this  case,  the  rate  of  diffusion  is  inde- 
pendent of  B.   Since  these  points  have  been  treated  ex- 

tensively  by  Grad  -^    ,   we  shall  not  elaborate  on  them 
__  — _ 

•^Harold  Grad,  Notes  on  Magneto -Hydrodynamics,   NY0-6[|.86, 
Nos.  IV  and   VII. 
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further.   Instead,  the  author  takes  pleasure  In  thank- 
ing Professor  Grad  for  many  useful  discussions  on  these 
topics. 

Basic  assumptions  and  equations.   We  consider  a  steady- 
state  flow  in  which 


I 


B 

>  X      all  quantities  de- 


->  u 


pend  on  x  alone,   x 
being  the  direction  of 
flow  perpendicular  to  B, 


Pigiore  1. 

We  assvime  that  the  following  conditions  prevail: 

1)  Departures  from  neutrality  are  so  small  that 
the  electrical  forces  acting  on  a  vol\ime  element  of  the 
fluid  may  be  neglected  as  compared  to  the  magnetic 
forces. 

2)  The  flow  is  isothermal,  hence,   p  =  ap  where 
p   is  the  pressure,   p,   the  density  and  a,   a  para- 
meter.  A  more  realistic  assvimption  would  be  that  ener- 
gy flow  takes  place.   However,  this  would  not  change 
the  qualitative  nature  of  the  results  although  quanti- 
tatively they  would  differ.   Isothermal  flow  may  even  be 
a  good  approximation  in  some  cases,  e.g.,  when  radiation 
effects  are  large. 

The  equations  describing  the  motion  of  the  fluid 
are  the  conservation  equations  of  mass  and  momentum  to- 
gether with  Maxwell's  equations  and  Ohm's  law. 
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Conservation  of  mass;   pu  =  M  or  pu  =  oM  (1) 

where   M  is  a  constant 
of  Integration. 

Conservation  of  momentum:   pu.^^u  +  grad  p  =  JxB  (2) 

Maxwell's  equations:   curl  E  =  0  (3) 

curl   B  =  p.J  ill) 

Ohm's  law:  E'  =  RJ  +  ^J  x  B  (5) 

where   E'  =  E  +  u  x  B 

Isothermal  flow:   p  =  ap  (6) 

We  note  that  R  =  -r  ,   hence,  is  independent  of  B 
whereas  E,     may  depend  on  B.   The  second  term  on  the 
right-hand  side  of  equation   (5)   is  the  Hall  vector 
term.   We  shall  now  show  that,  for  the  particular  ge- 
ometry considered  here,  this  term  drops  out  of  the  equa- 
tion of  motion  for  B.   The  latter  is  obtained  from 
equations   (3),  (i|)   and   (5).   From  c\irl   E  =  0  we 
obtain   -curl(uxB)  +  R  curl  J  +  C  curl  (JxB)  =  0        (7) 
But  curl (JxB)  =  0   since   J  only  has  a  y   component. 
Using  equations   (i^.)   and   (?)   and  div  B  =  0,   we 
obtain 

curl  curl  B  =  li  curl  J  =  •=  curl  (uxB)  (8) 

This  is  exactly  the  relation  obtained  when  Ohm's  law 
is  written  in  the  form  J  =  d(E  +  uxB),  without  the 
Hall  vector. 
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Equation  (8)  may  be  integrated,  yielding 


Bu  =  ^^  +  e  (9) 


where   e   Is  a  constant  of  Integration  equal  to  the  y 
component  of  the  electric  field.   By  using  equation 
([j.),   the  momentum  equation  may  also  be  expressed  In 
integrated  form: 

>2 


R 
Mu  +  p  +  -^  =  TT 

where   tt   is  another  integration  constant.   Equations 
(1),  (9)   and   (10)   constitute  the  basic  system  for 
determining  the  quantities  u,  p  and  B.   It  can  be 
used  to  describe  either  shock  flows  or  diffusive  phe- 
nomena.  The  difference  between  the  two  cases  shows  up 
in  the  number  of  parameters  that  are  required  in  order 
to  specify  the  problem.   This  may  be  understood  more 
easily  by  using  a  cylindrical  geometry  in  which  the 
magnetic  field  has  the  direction  of  the  axis  and  the 
flow  lies  in  planes  perpendicular  to  the  axis.   All 
variables  depend  on  r  alone,  namely,  the  distance 
along  the  radius  from  the  axis  of  the  cylinder.   The 
equation  tt  =  -  curl  E,   written  in  Integral  form  is 

\  4|  .  dS  =  -  1  curl  E  .  dS  =  -  ^  E  .  d/ 
A  A 
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where   d/^  is  an  element  of  length  along  the  curve  en- 
closing the  area  A,  which  is  perpendicular  to  and  con- 
tains the  axis  of  the  cylinder.   If  no  sources  are  enclosed, 

curl  B 
then  the  expression  on  the  right  vanishes  or  — — —  + 

+  uxB  =  grad  ^  .   Now,  grad  <-f>     lies  in  the   0  di- 
rection and  depends  on  r  alone.   In  order  to  satisfy 
these  requirements,  ^p  must  be  directly  proportional 
to  ©.   However,  the  integral  would  then  no  longer 
vanish.   The  only  way  to  satisfy  all  requirements  is 
for  ^  to  vanish  everywhere  hence  for  grad  C(?    to 
vanish.   In  the  case  of  shock  flows,  sources  are  needed 
along  the  axis  of  the  cylinder  in  order  to  maintain  a 
steady  flow  of  fluid  through  the  shock  front,  hence 
grad  Q^  =  e  7^  0,   In  the  case  of  diffusion,  however, 
this  is  not  necessary,  hence   grad  f  =  e  =  0.   For 
shocks,  therefore,  an  extra  parameter  (hence  boundary 
condition)  is  needed  in  order  to  specify  the  problem. 
This  case  will  be  treated  in  a  paper  to  follow  and  will 
not  be  pursued  here.   The  same  argument  can  be  extended 
somewhat  artificially  to  the  plane  case  by  allowing 
r  -->  00  and  the  sources  to  recede  to  infinity. 

Returning  to  diffusion,  it  would  appear  that 
three  parameters  are  necessary  in  the  problem,  namely, 
M,  TT  and  a.   However,  if  equations   (1),  (9)   and  (10) 
are  expressed  in  dlmensionless  form,  only  one  essential 
parameter  P  remains.   Introducing  the  dlmensionless 
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quantities      p  =  -^    ,      B  =  -^   ,      u  =  7-.  ,      x  =  /a"  ^Cm   , 
j5  =  Aj-—   ,      these   equations  become 


pu  =  ^  (11) 

^u  +  p  +  |-  =  1  (12) 

Bu  =  -^  (13) 

die 

ConclusioiB 

It  is  now  a  simple  matter  to  prove  the  assertions 
made  in  the  Introduction. 

1)  The  only  strict  scaling  of  a  flow  is  one  In 
which  P  is  maintained  constant.   In  such  a  case  we 
might  wish  to  compare  experiments  having  large  and  small 
magnetic  fields.   For  large  fields,   tt  would  have  to  be 
increased  in  order  to  keep  5  constant  hence  p  would 
also  have  to  be  increased  in  order  to  keep  p  constant. 
However  u  would  be  unaffected.   We  conclude  that   u 

is  independent  of   B   in  this  strict  sense. 

2)  Other  types  of  scaling  are  possible.   These, 
however,  are  only  local  and  allow  the  dependence  on  x 
to  change.   For  example   15  and  tr  might  be  allowed  to 
vary  keeping  -t^  constant.   It  is  then  easily  shown  that 

uOC— 3  for  large   B.   This  agrees  exactly  with  Chapman 

^  (3) 

and  Cowling's  result    which  is  derived  on  the  basis  of 

kinetic  theory  using  the  same  type  of  scaling. 


Chapman  and  Cowling.   Ibid,  pp.  327-328 • 
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3)      The   basic    system  of  equations  yields 


u  = 


-(1-  f-)   l/il'  1^)^-   k^^ 
15 


Hence,  for  a  given  flow,   u  is  in  general  not  propor- 
tional to  any  specific  power  of  B. 

I|.)   pu,   the  mass  flow,  is  independent  of  B 
since   B   depends  on  x  and   pu  does  not. 
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